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For the preparation of the experimental search for charged Higgs bosons at the LHC detailed stud- 
ies have been made to determine the expected charged Higgs boson production cross-sections and 
decay branching ratios at = 14 TeV. In the mass regime below the f -quark mass the expected 
production cross-sections are discussed using PYTHIA and FeynHiggs program packages based 
on the decay t — > H + b. For higher masses Next-to-Leading-Order (NLO) calculations have been 
used, and particular attention has been given to the intermediate-mass region. The decay branch- 
ing ratios have been studied with the program packages FeynHiggs and HDecay. Higher-order 
corrections (Aj, corrections) in the MSSM are consistently taken into account. Two benchmark 
scenarios are considered, one of them the 'mhmax'. 
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1. Introduction 

Charged Higgs bosons are naturally predicted in non-minimal Higgs scenarios, for instance in 
Two Higgs Doublet Models (THDMs), and specifically in the Minimal Supersymmetric Standard 
Model (MSSM) [1]. 

At present, a lower bound on the charged Higgs boson mass of about 80 GeV exists from 
LEP [2, 3], remarkably close to the previously simulated sensitivity [4]. With initial Tevatron 
data, upper limits were placed on BR(? — > H + b) for different charged Higgs decay scenarios [5, 
6]. Starting from 2009/2010, the LHC at CERN will enable the discovery or the ruling out of 
the existence of such a particle over a large portion of both the THDM and MSSM parameter 
space (masses up to around 400 GeV). The exact value of the reach depends on the value of tan/3 
(reviews [7-9] and a recent study [10]). The expected sensitivities for the LHC have been reported 
at this workshop [11, 12]. 

This comparative study focuses on the latest developments in the production cross-section and 
branching ratio determination, and is structured as follows. First the parameter points in the MSSM 
are defined in section 2. Then in section 3 the determination of the branching ratio BR(? — > H + b) 
is addressed and results from PYTHIA and FeynHiggs calculations are discussed. This branching 
ratio is used to determine the H + production cross-section in the low-mass region (section 4). 
Section 5 addresses the intermediate-mass region around m H + = 170 GeV. In the high-mass region 
the process gb — > H + is dominant and the calculations are described in section 6. The H + branching 
ratios are discussed in section 7. Systematic uncertainties are the focus in section 8. Section 9 
describes the data-base structure for storing the cross-section and branching ratio values. 

2. MSSM scenarios 

Two scenarios in the MSSM are considered. They are described by the following parameters. 
Scenario A: m, = 175 GeV, M SU sy = 500 GeV, A, = 1000 GeV, = 200 GeV, M 2 = 1000 GeV, 
M 3 = 1000 GeV. Scenario B ("mhmax") [13]: m t = 175 GeV, M SU sy = 1000 GeV, X, = 2000 GeV, 
where A, = X, + pi / tan /3 , = 200 GeV, M 2 = 200 GeV, M 3 = 800 GeV. 

The A/, corrections are calculated in FeynHiggs v2.6.2 [14] for these two cases in the H + 
couplings. For tan/3 = 50 they are A/, = 0.45 for scenario A, and Ab = 0.21 for scenario B. The Ab 
corrections modify the ft-quark mass m ° orrected =mb /( \ -\-Ab) [15]. 

3. t — >• H+b branching ratios 

The BR(? — > H + b) values have been determined with PYTHIA v6.325 [16] and FeynHiggs 
v2.6.2 [14] . In FeynHiggs the formula from Ref. [17] is implemented, and furthermore it includes 
the Ab corrections depending on the MSSM parameters [15]. The computations have been per- 
formed for MSSM scenarios A and B. An example is shown in Fig. [I] for scenarios A as a function 
of tan /3 . The FeynHiggs calculations include A^ corrections, while the PYTHIA calculation does 

'in this study the branching ratios were also produced with FeynHiggs v2.6, however, then discarded as the differ- 
ences in version 2.6 and 2.6.2 were only attributed to a programming correction ('bug fix') in the latter version. 



2 



Cross-sections and branching ratios for charged Higgs searches 



Andre Sopczak 



not. The At, corrections are positive and grow with tan j8 which explains the observed increase in 
BR(? — > H + b) difference with increasing tan/3 (Fig. [I]). 

The first step in the determination of the production cross-section is the calculation of the 
branching ratio BR(? — > H + b) in the low mass region 90 to 170 GeV. These branching ratios are 
also shown in Fig. [I] for scenario A, calculated with FeynHiggs. 




Figure 1: Left: Expected branching ratio BR(f — > H + b) for MSSM scenario A. Right: Expected branching 
ratio BR(f — > H + b) for MSSM scenario A as described in the text, calculated with FeynHiggs version 2.6.2. 
For very large tan j3 values and low charged Higgs masses (90 GeV), the model enters a non-perturbative 
regime and loop calculations are no longer valid, as indicated by the ending of the 90 GeV curve. 



4. Low-mass cross-section 

The charged Higgs boson production cross-section has been evaluated with different program 
packages for low-mass and high-mass regions. If the charged Higgs boson mass m H + satisfies 
m H + < m t — nib, where nib is the ^-quark mass, H + particles could be produced in the on-shell 
process t — ► bH + , the latter being in turn produced from gg fusion and qq annihilation. Contri- 
butions of gluon fusion to the production of charged Higgs at hadron colliders were pointed out 
previously [18]. This approximation has customarily been used in event generators for m H + < m t . 

Charged Higgs production is denoted by qq, gg — > ft — > tbH + if due to (anti-) top decays and 
by qq, gg — > tbH + if further production diagrams are included [19]. Owing to the large ?-quark 
decay width (T t ~ 1.5 GeV) and to the additional diagrams which do not proceed via direct ft 
production [20-22], charged Higgs bosons could also be produced at and beyond the mass threshold 
(m t — rtib) in the tbH + process. The importance of these effects in the so-called 'threshold' or 
'transition' region (intermediate-mass region m H + m t ) was emphasized in various Les Houches 
proceedings [23,24] as well as in Refs. [19,25-27], and the calculations of Refs. [20,21] (based 
on the appropriate qq,gg —> tbH + description) are now implemented in HERWIG [28-31] and 
PYTHIA [16, 19]. A comparison between the two generators was carried out in Ref. [19]. In 
addition, in the mass region near the f-quark mass, a matching of the calculations for the qq, gg — > 
tbH + and gb — ► tH + processes is required already in leading order [32]. In the kinematic region 
where t — ► H + b is possible, this process could dominate the charged Higgs production. 

The cross-sections in the low-mass region for the charged Higgs boson masses 90, 110, 120, 
130 and 150 GeV have been calculated from the higher-order improved ft cross-section {o t j = 
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833 pb [33]) and the BR(? -> H + b) determined from FeynHiggs (version 2.6.2) [14]: a tbH + = 
2 ■ a ff BR(t -> H + b) [1 - BR(? -> H + b)} . Results are shown in Fig. | (right plot). 

5. Intermediate-mass cross-section 

Charged Higgs bosons with a mass of 170 GeV would predominantly be produced by the gb — > 
tH + process. The intermediate-mass region has been studied in NLO [34]. The f-quark mass has 
been fixed to 175 GeV. For this scenario with a 5 GeV mass difference between charged Higgs and 
?-quark masses, the additional cross-section from the t — ► H + b process amounts to an increase of 
about 20 to 30%. We have taken into account this increase in the derivation of the production cross- 
section by adding both cross-sections. The cross-section increase depends strongly on the mass 
difference between charged Higgs and ?-quark masses, and also on the treatment of the running 
fe-quark mass. Results are shown in Fig. || (right plot, 170 GeV curve). 

6. High-mass cross-section 

At hadron colliders, the main contribution to charged Higgs boson production is through the 
twin processes gg — ► tbH + and gb — > tH + for m H + > m t . These are called twin processes since they 
correspond to two different approximations describing the same basic process. For charged Higgs 
boson masses above the ?-quark mass, the 2^2 process is dominant, due to the resummation of 
potentially large logarithms in the &-quark parton density [35]. In the high-mass region the Next- 
to-Leading Order (NLO) production cross-section calculation is applied [34,35]. In this case, the 
parton shower produces an outgoing ft-quark of relative small transverse momentum. In the region 
of phase space where the outgoing ft-quark has large transverse momentum, the parton shower does 
not give a good description of the process, and the full 2^3 description is needed. However, these 
two descriptions overlap for small transverse momenta of the ft-quark, necessitating a matching 
procedure to remove resulting double counting [32]. For charged Higgs boson masses below the 
?-quark mass, the 2^3 process dominates since it incorporates on-shell ?-quark pair production 
with subsequent decay into a charged Higgs boson. Matchig [36] is a new leading-order event 
generator based on the work presented in Ref. [32] which matches the two processes by producing 
negative weight events from an identified double-counting term. In the high-mass region the NLO 
program package [34] also avoids double counting. 

For the charged Higgs boson masses 200, 250, 350, 400, 500 and 600 GeV, the BR(f -> 
H + b) is kinematically suppressed and the cross-sections have been determined in NLO from the 
gb — > tH + process alone [35]. We have explicitly not calculated the dependence of the cross- 
sections on the MSSM parameters with the NLO program package [34, 35] 2 . These higher-order 
corrections can be large, depending on the MSSM scenario. They depend primarily on tan/3. In 
order to determine these corrections, first the Ab corrections are calculated with the FeynHiggs 
package [14]. Then a reduction factor / = 1/(1 + Ab) 2 is calculated for the production cross- 
section. The previously determined NLO cross-sections are multiplied by this reduction factor. 

2 Recently, the charged Higgs production process at NLO has been implemented in the program package 
PROSPINO 2.1 including A b corrections. 
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The reduction factors for scenario A are shown in Fig. § as a function of tan/3. The figure shows 
also the cross-sections after the application of the reduction factors 3 . 
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Figure 2: Left: Reduction factors / which have been applied to the NLO charged Higgs boson production 
cross-section in the MSSM for scenario A. The reduction factors result from Aj, corrections and have been 
calculated with FeynHiggs v2.6.2 [14]. Right: Expected charged Higgs boson production cross-section in 
the MSSM for scenario A. For charged Higgs boson masses of 170 GeV and above the reduction factors / 
have been applied to the NLO cross-section calculations [35]. This allows a consistent treatment of the 
corrections for the production cross-sections and branching ratios over the whole mass range. 



7. H + branching ratios 

This section describes the H + branching ratio studies, in particular focusing on the off-mass- 
shell effects and the Ay corrections. The BR(H + —>...) values have been determined with the 
FeynHiggs v2.5, v2.6.2 and HDecay 3.2 program packages. FeynHiggs v2.5 does not include off- 
mass-shell effects while FeynHiggs v2.6.2 does. The difference between FeynHiggs v2.6.2 and 
HDecay 3.2 4 is that FeynHiggs includes the Ab corrections. 

While for m H + < m t charged Higgs bosons decay predominantly into a T-lepton and a neu- 
trino, or into a cs-quark pair; for m H + > m t both H + — ► tv t and H + — > tb are important decay 
channels. In the experimental search, the latter is much harder to disentangle from background 
than the former. The associated ?-quark decays predominantly into a W boson and a ft-quark. 

The branching ratios have been determined with the FeynHiggs [14] and HDecay [38] pro- 
gram packages. The detailed comparison showed very good agreement between FeynHiggs and 
HDecay calculations for the branching ratios in the low-mass region (m H + < m t ), including the 
virtual effects which lead to H + — > tb contribution in this mass region 5 . In the high-mass region, 
vertex corrections (A/, terms) which are included in FeynHiggs lead to a significant variation of the 
branching ratio with tan j3 while the branching ratios calculated with HDecay are largely indepen- 
dent of tan j8 in the high-mass region. The branching ratio comparison is given in Table [I]. 

Examples of branching ratios are shown in Fig. |3] in the low-mass region for a 130 GeV and in 
the high-mass region for a 600 GeV charged Higgs boson [14]. The dependence of the branching 

3 However, for any sub-dominant decay channel of a heavy charged Higgs boson, such as H + — > T + v, Aj, corrections 
cancel to a large extent [37]. 

4 Since version 3.3 Aj, corrections are included in HDecay for the charge Higgs sector. 

5 These virtual effects were only included in FeynHiggs version 2.6, release July 2007. In previous versions of 
FeynHiggs the tb branching fraction due to virtual effects was zero. 
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170 
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400 


tan/3 = 3, FH 


0.901/0.766 


0.009/0.008 


0/0.1485 


0.004/0.004 


0.000/0.000 


0.978/0.978 


3, HD 


0.745 


0.008 


0.133 


0.004 


0.000 


0.984 


10, FH 


0.990/0.988 


0.006/0.006 


0/0.002 


0.146/0.146 


0.001/0.001 


0.845/0.845 


10, HD 


0.974 


0.006 


0.001 


0.112 





0.88 


60, FH 


0.991/0.991 


0.006/0.006 


0/0.000 


0.336/0.336 


0.002/0.001 


0.660/0.662 


60, HD 


0.976 


0.006 





0.143 





0.854 



Table 1: Branching ratios BR(H + — > TV, cs, tb) for FeynHiggs (FH), left: v2.5, right v2.6.2, and HDecay 
(HD) v3.2. Version 2.6.2 includes off-mass-shell effects, while version v2.5 does not. It has been explicitly 
checked that the difference between FeynHiggs v2.6.2 and HDecay are due to the Aj corrections [39]. 



ratio as a function of the mass is also shown in Fig. [|. When kinematically allowed the decay 
H + — > X°X + could become large, as calculated with FeynHiggs v2.6.2. This decay mode is not 
addressed in the current discovery analyses and this branching fraction is not discussed further. 




tanp H* mass (GeV) 



Figure 3: Left: Expected charged Higgs boson branching ratios in the MSSM for scenario A for an example 
of a light and a heavy charged Higgs boson [14]. Right: Expected charged Higgs boson branching ratios in 
the MSSM for scenario A as a function of the charged Higgs boson mass [14]. 

8. Systematic Uncertainties 

The systematic uncertainties are discussed for the charged Higgs branching ratio BR(t — > H + b) 
and production cross-section in the low-mass region for m H + < m t . The charged Higgs boson 
branching ratios BR(// + — > xv,cs,tb) have been determined with the same program package and 
similar systematic uncertainties apply. The production cross-sections in the high-mass region have 
been determined with a NLO program package [35], and the systematic uncertainties are discussed 
separately. 

The values for BR(f — > H + b) and BR(// + — > %v,cs,tb) have been calculated with FeynHiggs 
v2.6.2 which includes vertex corrections in the framework of the MSSM (Aj terms). Systematic 
uncertainties from higher-order loop corrections to the tbH + vertex and due to the running of the 
c and s-quark masses are expected [39]. Upper limits on the resulting systematic uncertainties are 
conservatively estimated to be: 

AT(f -» H+b)/F < 10%, AT(H+ -» tv)/T < 5%, AT(H+ -> cs,tb)/T < 10%. 
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Systematic uncertainties on the charged Higgs boson production cross-section in the high- 
mass region can occur primarily from two sources. First from the renormalization scale and fac- 
torization scale dependence, and second from the fact that supersymmetry loop corrections are not 
(yet) included in the calculations [40,41]. The one-loop contributions largely improve the the- 
oretical uncertainty of the leading order (LO) cross-section. The remaining uncertainty can be 
estimated from the scale dependence. The variation of the gb — ► tH + production cross-section for 
0.1 < jU/jU C entrai < 10 was considered, where /i C entrai = ( m t + Wff+)/5. The resulting systematic 
uncertainty on the production cross-section is below 20%. 

Corrections from supersymmetric particles in the MSSM are not included in the NLO produc- 
tion cross-section calculations [40,41]. In this study, these supersymmetry loop corrections (Af, 
corrections) have been taken into account independent of the NLO calculations using the Feyn- 
Higgs package and thus no additional uncertainty beyond the 20% is assigned 6 . 

9. Database 

The branching ratio BR(? — ► H + b), the charged Higgs boson production cross-section in the 
low-mass, intermediate-mass and high-mass regions, the A^ values and corresponding cross-section 
reduction factors, and the H + branching ratios have been determined for various charged Higgs 
boson masses between 90 and 600 GeV and for tan j8 values between 1 and 90. The calculated 
values have been stored in a database in root format for scenarios A and B. 

10. Conclusions 

Comparative studies for the expected charged Higgs boson production cross-sections and 
branching rations have been performed for searches in the initial LHC data. The production cross- 
sections in the low-mass region and all charged Higgs branching ratios have been calculated with 
FeynHiggs v2.6.2. The production cross-sections in the high-mass region have been determined 
with a dedicated NLO program package and the dependence on MSSM parameters has been added 
using the FeynHiggs package. 
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